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Abstract 



These notes are a summary of lectures given at the Instituto de Astronomia of 
the Universidad Nacional Autonoma de Mexico (UN AM), the Dipartimento di 
Fisica of the Universitd degli studi di Salerno (Italy), and the Instituto de Fisica 
Corpuscular of the Universitat de Valencia (Spain) during the year 2012. 
Basic mathematical aspects of Palatini theories of gravity, which are constructed 
assuming that metric and connection are independent geometrical entities, are 
briefly introduced and discussed, and followed by a detailed derivation of the 
field equations of some rather general Lagrangian theories. Applications to the 
early universe are considered paying special attention to the avoidance of the 
big bang singularity. The analysis of non-singular cosmologies carried out in 
arXiv:1005.4136 [gr-qc] is extended to include several perfect fluids. 



Contents 



1.1 Introduction 1 

1.2 Diffcrcntiablc manifolds, affine connections, and the metric. ... 3 

1.3 Dynamics of Palatini theories 6 

1.3.1 Example: f(R,Q) theories 8 

1.3.2 Volume-invariant and torsionless f(R, Q) 

1.4 f(R, Q) theories with a perfect fluid 

1.4.1 Workable models: f(R, Q) = f(R) + aQ 

1.5 Nonsingular cosmologies in f(R,Q) theories 

1.5.1 Homogeneous cosmologies in f(R,Q) theories 

1.5.2 Shear 

1.5.3 Expansion 

1.5.4 Limit to f(R) 

1.5.5 Bouncing f(R) Cosmologies 

1.5.6 Nonsingular Universes in f(R, Q) Palatini theories .... 

1.5.7 Anisotropic Universe 

1.5.8 Example: radiation universe 

1.6 Conclusions and open questions 



1.1 Introduction 

The impact of Einstein's fundamental idea of gravitation as a curved space-time 
phenomenon on our current understanding of the Universe has been enormously 
successful. A key aspect of his celebrated theory of General Relativity (GR) is 
that the spatial sections of four dimensional space-time need not be Euclidean. 
The Minkowskian description is just an approximation valid on (relatively) local 
portions of space-time. On larger scales, however, one must consider deforma- 
tions induced by the matter on the geometry, which must be dictated by some 
set of field equations. In this respect, the predictions of GR are in agreement 
with experiments in scales that range from millimeters to astronomical units, 
scales in which weak and strong field phenomena can be observed [ |Will, 1993| ]. 
The theory is so successful in those regimes and scales that it is generally ac- 
cepted that it should work also at larger and shorter scales, and at weaker and 
stronger regimes. The validity of these assumptions, obviously, is not guaran- 
teed a priori regardless of how beautiful and elegant the theory might appear. 
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Therefore, not only must we keep confronting the predictions of the theory with 
experiments and/or observations at new scales, but also we have to demand 
theoretical consistency with the other physical interactions and, in particular, 
in the quantum regime. 

For the above reasons, we believe that scrutinizing the implicit assumptions 
and mathematical structures behind the classical formulation of GR could help 
better understand the starting point of some current approaches that go be- 
yond our standard model of gravitational physics. At the same time, this could 
provide new insights useful to address from a different perspective some current 
open questions, such as the existence of black hole and big bang singularities 
or the cosmic speedup problem. In this sense, Einstein himself stated that 
" the question whether the structure of [the spacetime] continuum is Euclidean, 
or in accordance with Riemann's general scheme, or otherwise, is ... a phys- 
ical question which must be answered by experience, and not a question of a 
mere convention to be selected on practical grounds" [ |Feigl fc Brodbeck, 1953] ]. 
From these words it follows that questioning the regime of applicability of the 
Riemannian nature of the geometry associated with the gravitational held and 
considering more general frameworks are legitimate questions that should be 
explored by all available means (theoretical and experimental) . These are some 
of the basic points to be addressed in this work. 

In this chapter we explore in some detail the implications of relaxing the 
Riemannian condition on the geometry by allowing the connection to be deter- 
mined from first principles, not by choice or convention. This approach, known 
as metric-afhne or Palatini formalism [ |01mo, 201 la| ] , assumes that metric and 
connection are equally fundamental and independent geometrical entities. In 
consequence, any geometrical theory of gravity formulated in this approach 
must provide enough equations to determine the form of the metric and the 
connection (within the unavoidable indeterminacy imposed by the underlying 
gauge freedom). We derive and discuss the field equations of a rather general 
family of Palatini theories and then focus on two particular subfamilies which 
have attracted special attention in recent years, namely, f(R) and f(R, Q) the- 
ories. The interest in studying these particular theories lies in their ability to 
avoid (or soften in some cases) big bang and black hole singularities and their 
relation with recent approaches to quantum gravity. Here we will focus on the 
early-time cosmology of such theories. 

The content is organized as follows. We begin by briefly reviewing in section 
11.21 the basics of differentiable manifolds with afhne and metric structures, to 
emphasize that metric and connection are equally fundamental and indepen- 
dent geometrical objects. In section 11.31 a derivation of the field equations for 
a generic action depending on the metric and the Riemann tensor is presented 
taking into account also the presence of torsion. In section 11.41 we discuss a 
particular family of Lagrangians of the form f(R, R^ V R^ U ) in combination with 
perfect fluid matter, and prepare the notation and field equations needed to 
study the dynamics of those theories. We then focus on the early-time char- 
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acteristics of isotropic and anisotropic homogeneous cosmologies 11.51 and show 
that nonsingular bouncing solutions exist for f(R) and f(R, Q) models (subsec- 
tions [TITS] and [T3TS1 respectively). We conclude with a discussion of the results 
presented and point out some open questions that should be addressed in the 
future. 



1.2 Different iable manifolds, affine connections, 
and the metric. 

In this section we quickly review some of the mathematical structures needed 
to construct a geometric theory of the gravitational interactions. The goal is 
to put forward that metric and connection are equally fundamental and in- 
dependent geometrical entities, an aspect usually overlooked in the construc- 
tion of phenomenological extensions of GR. We will thus be more sketchy 
than mathematically accurate. For a more exhaustive and precise discussion 
of these topics see your favorite book on differentiable manifolds (or, for in- 
stance, [[Naka hara, 1990| ]). 



In the geometric description of gravitational theories, one begins by identi- 
fying physical events with points on an n-dimensional manifold M. The next 
natural step is to provide this manifold with a differentiable structure. One then 
labels the points p € M with a set of charts (Ui,<pi), where the Ui are subsets 
of M and ipi are maps from Ui to K™ (or an open subset of K n ) such that every 
p G M lies in at least one of the charts (Ui,<fi). If for any two charts (Ui,tpi) 
and {Uj,ipj) that overlap at some nonzero subset of points the map ipi o ipj 
is not just continuous but differentiable, then we say that M. is a differentiable 
manifold. 

Since the Euclidean view of vectors as arrows connecting two points of the man- 
ifold is not valid in general, to get a consistent definition we need to introduce 
first the concept of curve and tangent vector to a curve at a point. We thus say 
that a smooth curve j(t) in M. is a differentiable map that to each point of a seg- 
ment associates a point in M, 7(<) : t G [0, 1] — > M. In a chart (U, <p), the points 
of the curve have the following coordinate representation: x — <*p(jpt) — <f ° lif)- 
If we consider now a function /onM, where / is a map that to every p € A4 
assigns a real number (/ : M. — > M), the rate of change of / along the curve j(t) 
using the coordinates of the chart (U, tp) is given by 

ff(y°7(*)) _ M _ df dx»(t) _ df 

dt dt dx» dt [ ) dx» ' 1 ' 

where we have defined the components of the tangent vector to the curve in 
this chart as X M = dx^(t)/dt. Vectors can thus be seen as differential operators 
X = X^dfj, whose action on functions is of the form X[f] — X^d^f, thus provid- 
ing a natural notion of directional derivative for functions. The set {e M = d^} 
defines a (coordinate) basis of the tangent space of vectors at the point p, which 
we denote T p A4. Obviously, vectors exist without specifying the coordinates. 
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Under changes of coordinates, we have V = V^e^ = V a e a — U^lfs-e^, which 
implies the well-known transformation law V 1 = V a ^^ for the vector compo- 
nents. 

When a vector is assigned smoothly to each point of A4, it is called a vector 
field over A4 . Each component of a vector field is thus a smooth function from 
M — > M. Given a vector field X, an integral curve of X is defined as the 
curve whose tangent vector coincides with X. For infinitesimal displacements 
of magnitude e in the direction of X, a given point p of coordinates becomes 
<t^(x) = x^ + eX^(x). This transformation also induces a correspondence be- 
tween vectors of the tangent spaces T X M. and T a r x )M. The effect of these 
transformations on a vector field Y(x) leads to the concept of Lie derivative, 
whose action on vector fields is defined as 

C X Y = [X v d v Y»{x) - Y v d v X»{x)] e M = [X, Y] . (1.2) 

This derivative operator is independent of the choice of coordinates and follows 
naturally from the differential structure of the manifold. It satisfies a number of 
useful properties such as bilinearity in its two arguments, C X (Y + Z) = CxY + 
C X Z, Cx+yZ = C X Z + C Y Z, and the chain rule C x fY = (C x f)Y + f£ x Y, 
with Cxf = X[f\. Though the Lie derivative provides a natural directional 
derivative for functions, it docs not work in the same way for vectors and ten- 
sors of higher rank. In fact, since the partial derivatives of the vector X appear 
explicitly in CxY , two vectors whose components at a given point have the same 
values but whose partial derivatives at the point differ do not yield a vector that 
points in the same direction, i.e., they are not proportional. Therefore, in order 
to introduce a proper notion of directional derivative for vectors and tensors, 
we need to introduce a new structure called connection which specifies how 
vectors (and tensors in general) are transported along a curve. 

Manifolds with a connection. We are thus going to introduce a derivative 
operator, which we denote by V, such that given two vector fields X and Y we 
obtain a new vector field Z defined by Z = V X Y. This derivative operator 
must be bilinear in its two arguments, V x(Y + Z) = V xY + Vj2, V x+yZ = 
V X Z + V Y Z, must satisfy the chain rule V x (fY) = {V x f)Y + fV x Y, with 
Vxf = X[f], and must also behave as a natural directional derivative in the 
sense that VfxY — /V xY to guarantee that any two proportional vectors 
yield a result that points in the same direction. In a given coordinate basis, 
we have V X Y = X^V ^{Y" e v ) = X^ {e^\Y v ]e v + Y u V efi e v ). If our manifold 
is to— dimensional, defining m 3 functions called connection coefficients by 
V efi e v = T*„e\ we find that the last requirement, V/jfF = f\7 x Y, is naturally 
satisfied. We thus find that 



V X Y = X^ 



A 



— + r A y v 

dx» ^ 



ex • (1.3) 



The connection coefficients specify how the basis vectors change from point to 
point and, in principle, can be arbitrarily defined. Under changes of coordinates, 



4 



these coefficients transform as follows: 



which implies 



dx a 





dx a 




~ dxf 1 



dx x <9 2 x 7 
dx v dx x dx v 



UJj p7 
Q x u af. 



\ dx x dx a dxP~~ dx x d 2 x^ 

tl = — — ^-n fl + . 



(1.4) 



(1.5) 



M!y ftr 7 da^ ftr* 9iT dx»dx v ' 

This transformation law indicates that the connection coefficients do not trans- 
form as tensorial quantities. Therefore, the connection cannot have an intrinsic 
geometrical meaning as a measure of how much a manifold is curved. As intrin- 
sic geometric objects, we can define the torsion tensor 



T(X,Y)=V X Y -V Y X -[X,Y] , 
and the Riemann curvature tensor 

R(X,Y)Z = \7 X \7 Y Z- \7 Y \7 X Z- V [X , Y ]Z ■ 
In a coordinate basis, these tensors have the following components: 



-R(e M ,e„)eA = 



d T p 



d 



(1.6) 

(1.7) 

(1.8) 
(1.9) 



With the introduction of the connection, one can define the notion of parallel 
transport. Given a curve j(t) such that its tangent vector in a given chart has 
coordinates X^ = dx* 1 (t) /dt,we say that a vector Y is parallel transported along 
7(i) if V X Y = 0. In components, this equation reads ^f- + dx d ^ Y 13 = 0, 
where d/dt = X^d^. Geodesies are defined as those curves which are parallel 



or 



r M X a X P = 0. 



transported along themselves, namely, V X X — , -.. 

Manifolds with a metric. So far we have been able to construct a number 
of geometrical objects such as directional derivatives of vectors and tensors in 
general, the torsion and Riemann tensors, geodesic curves, . . . without the need 
to introduce a metric tensor. A metric tensor provides a notion of distance 
between nearby points and allows, among other things, to determine lengths, 
angles, areas, and volumes of objects which are locally defined in space-time. 
Formally, a (pseudo-Riemannian) metric tensor is a symmetric bilinear form that 
at each p e M. satisfies g p (U, V) — g p {V, U) for any two vectors U, V e T P M. 
and g p (U, V) = for any U € T p M iff V = 0. The metric tensor allows to define 
an inner product between vectors and also gives rise to an isomorphism between 
T p M and the dual space of one- forms T*M. In a coordinate basis, it can be 
represented by g = g ^ v dx^ ® dx v , where the differentials dx^ form a basis of 
T*M. In manifolds with a metric, one can impose a particular relation between 
the metric and the connection by demanding that the scalar product of any 
two vectors which are parallel transported along any curve remains covariantly 
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constant. This condition can be translated intc0 S7fj,g a /3 = 0, which implies that 
(recall that T*, = V 0a - I*,) 

2rf^) + (T^g^ + TP a g pv ) g° x = g x ? [d^ + d v g m - d p g^} . (1.10) 

From the right-hand side of this equation, one defines the Levi-Civita connection 
as 

L i±v = '— [d^gpu + d u g Pfl - dpg^} . (1.11) 

From this definition it follows that when the torsion vanishes, the connection is 
symmetric and coincides with the Levi-Civita connection. In that case, when 
= L x , we say that the associated geometry is Riemannian. It should be 
noted that though connections are not tensors, the difference between any two 
connections is a tensor. This, in particular, allowed us to construct the torsion 
tensor. With more generality, when the manifold is provided with a metric, any 
connection can be expressed as 

r^ = i^ + A^, (i.i2) 

where A x is a tensor (which needs not be symmetric in its lower indices). 
Therefore, Palatini theories of gravity, in which metric and connection are re- 
garded as independent fields, can be seen as theories in which an additional 
rank-three tensor field A x v has been added to the gravitational Lagrangian. 



1.3 Dynamics of Palatini theories 

From the above quick review of the properties of differentiable manifolds with 
metric and affine structures, it is clear that metric and connection are equally 
fundamental and independent geometrical entities. In the construction of the- 
ories of gravity based on geometry, we will thus assume this independence and 
will require those theories to yield equations that allow to determine both the 
metric and the connection and the possible relations between them. For sim- 
plicity, we will assume that the matter is only coupled to the metric (which is 
consistent with the experimental tests of the equivalence principle [ Will, 1993] ]) 



but will allow an independent connection to appear in the gravitational sector 
of the theory. As pointed out above, this is equivalent to having, besides the 
metric, a rank-three gravitational tensor field. From a geometric perspective, 
this possibility seems much more natural and fundamental than considering, for 
instance, scalar fields in the gravitational sector, though scalar-tensor theories 
have traditionally received much more attention in the literature. 

We begin by deriving the field equations of Palatini theories in a very gen- 
eral case and then consider some simplifications to make contact with the lit- 
erature. For a generic Palatini theory in which the connection appears through 
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the Riemann tensor or contractions of it, the action can be written as follows 
[ |Qlmo, 201lE| 1 



(1.13) 



where S m is the matter action, ip represents collectively the matter fields, k 2 is 
a constant with suitable dimensions (if f = R, then k 2 = 8irG), and 



R° 



r) V a 



1 fi\ L v& 



pa 

1 i/A 1 nfi 



(1.14) 



represents the components of the Riemann tensor, the field strength of the 
connection T^. Note that since the connection is determined dynamically, i.e., 
we assume independence between the metric and affine structures of the theory, 
we cannot assume any a priori symmetry in its lower indices. This means that 
in the variation of the action to obtain the field equations we must bear in mind 
that Tp 7^ r"o, i.e., we admit the possibility of nonvanishing torsion. It should 
be noted that in GR energy and momentum are the sources of curvature, while 
torsion is sourced by the spin of particles [[?]]. The fact that torsion is usually 
not considered in introductory courses on gravitation may be rooted in the 
educational tradition of this subject and the fact that the spin of particles was 
discovered many years after the original formulation of GR by Einstein. Another 
reason may be that the effects of torsion are very weak in general, except at very 
high densities, where the role of torsion becomes dominant and may even avoid 
the formation of singularities (see [ |Poplawski, 2010] ] for a recent discussion and 
earlier literature on the topic). For these reasons, and to motivate and facilitate 
the exploration of the effects of torsion in extensions of GR, our derivation of 
the field equations will be as general as possible (within reasonable limits). 
We will assume a symmetric metric tensor g^ u = g vll and the usual definitions 
for the Ricci tensor R^ v = i? p MP „ and the Ricci scalar R = g^R^. The 
variation of the action (|1.13[) with respect to the metric and the connection can 
be expressed as 



SS 



1 

2^2 



Straightforward manipulations show that SR a can be written as 

SBf Buu = V„ (STZp) - V, (5T^) + 2S*JT^ , (1.16) 



0HU 



where = (r* — T„ fl )/2 now represents the torsion tensor tensor (note the 
additional 1/2 factor as compared to our initial definition in ll.8[) . From now on 
we will use the notation P Q! ^" / = g^J^ — ■ In order to put the 5R a p^v term in 
(|1.15[) in suitable form, we need to note that 



d 4 x 



-gJ") 
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where = Pj^ST^. Since, in general, V^yf^gJ'*) = d^y/^gJ^) 
2S°„T l /—gJ fl , we find that (|1.17[) can be written as 

Ir=J d'x [d,(V^Jn - 5T« {v, (y=gPj»») 2S^V=gP a f »"} 

Using this result, (|1 . 15[) becomes 



(1.18) 



1 



'(|^-^) V- + ^U'-//-/") 



(1.19) 



We thus find that the field equations can be written as follows 
df f 



U0 _ 



K 2 H a 



1 



(1.20) 



-<J 



SZ p Pj° p + ZS^PjW (1.21) 



where PjW = (P a ^-P a ^)/2, = - 
tensor of the matter, and H a v ^ = 1 >>J 



2 SS„. 

1—9 $9>" 



is the energy-momentum 



; <5r° 



represents the coupling of matter 

to the connection. For simplicity, from now on we will assume that H a v ^ = 0. 
Eq. (|1.21[) can be put in a more convenient form if the connection is decom- 
posed into its symmetric and antisymmetric (torsion) parts, V^ v = C™ v + S™ v , 
such that V M A, = d„A v - C% v A a - S™ u A a = V°A V - S™ v A a and V^^/^ = 
Vj;^- S^aV-g- B y doin S tnis > fTTT^T]) turns into 



K 2 H a ^ = - 



(V^9PJ l ^ ] ) + S^PyfW _ S^P^W . (1.22) 



1.3.1 Example: f(R,Q) theories 

Eqs. (|1.20[) and (ll.22j) can be used to write the field equations for the met- 
ric and the connection for specific choices of the Lagrangian f(g f _ l ^,R a i3 llv )- 
To make contact with the literature [ |01mo, 2 011c ], [[Barraga n &: Olmo, 2010| ] 
Olmo et al., 2009| ], we now focus on the case f(R, Q) = fig^R^g^g^R^R^)- 



For this family of Lagrangians, we obtain 

Pj^ = S^MP" = S a » (f R g^ + 2f Q R^) , 



(1.23) 



where fx = dxf- Inserting this expression in (ll.22[) and tracing over a and u, 
we find that V^[y/^gM px ] = (2^/^/3)[S^M< 3x + (3/2)5? M A "]. Using this 



result, the connection equation can be put as follows 

1 -,C 



V- [V=gM^] = S v aX M^ - S% X M X » - S x aX M^ + -S'S^M^ (1.24) 
9 " 
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The symmetric and antisymmetric combinations of this equation lead, respec- 
tively, to 



1 



-9 



S v aX M^ - Sf^MW - S X X M^ + % (5^M px + 5^M vX \l.2h) 
= S u aX M {J3X) - S ( l x N& x) ~ S X X M [ ^ ] + ^ {5 u a M px - 5 a M uX ) (.1.26) 



Important simplifications can be achieved considering the new variables 



\JbV 



and taking the parameter a = 2/3, which implies that = f 



[ixv] 



(1.27) 
is such that 



S% v = 0. The symmetric and antisymmetric parts of T^ v are related to those 



ofT^ by 



•A 



Ci 



~{6ZS^ + 6 X SZ„) 

1 

A qct 



(1.28) 
(1.29) 

Using these variables, Eqs. (jl.25l) and (jl -26[) take the following compact form 



Sfiiv + 2 v San ^n^Zv 



1 



7C 



9^ 



Kp] 



/ p M, 



(rep) 



(1.30) 
(1.31) 



In these equations, = f R g pu + 2f Q R^ u \T), and = 2f Q R^(T), 

where R {M (T) = R (M (T) and R m (T) = R [M {Y) + § (fyS^, - 



In the recent literature on Palatini theories, only the torsionless case has 
been studied in detail. When torsion is considered in f(R) theories, Eqs. (|1.30j) 
and (|1.31[) recover the results presented in [ |Qlmo, 201 la| ]. In general, those 
equations put forward that when the traceless torsion tensor S^ x vanishes, the 
symmetric and antisymmetric parts of M$ v decouple. The dynamics of these 
theories, therefore, can be studied in different levels of complexity. The simplest 
case will be studied here and consists on setting S v aX and R\^} to zero. A more 
detailed discussion of the other cases can be found in [ |01mo fc Rubiera, 2012| ]. 

1.3.2 Volume-invariant and torsionless f(R,Q) 

When the torsion is set to zero, it can be shown [ |Schouten, 1951| ], [ |Hehl et al., 1 995 ] 
that the vanishing of i?[ M „] guarantees the existence of a volume element that 
is covariantly conserved by . The rank- two tensor that defines that volume 
element must be a solution of (ll.30|) . which in this case takes the form 

< [V=9 {fR9 f3u + 2/ Q i?^(r))] = . (1.32) 
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Note that here R^ V (T) is symmetric because we are taking R^„](T) = 0. To 
obtain the solution of (jl.32l) , we first consider (11.201) particularized to our theory 
(with S£ x and R\p, v ] set to zero), 

fRRfiiJ — + ^fQRfiaR" y = ^Tp,v , (1.33) 

and rewrite it in the following form 

f R B/ - Is/ + 2f Q B tl a B a v = K 2 T/ , (1.34) 

where we have defined B/ = i? MQ g ai/ . This equation can be seen as a second- 
order algebraic equation for the matrix B, whose components are [B]^ = B/ . 
The solutions to this equation imply that B is an algebraic function of the 
components of the stress-energy tensor T fl v ) i.e., B = B(T). This relation is 
very important because it allows to express f|l .32[) in the form 

< [V=99 PX (f R SX + 2f Q B x ")] = , (1.35) 

where now Jq and B a u are functions of the stress-energy tensor of the 
matter. The connection, therefore, can be obtained by elementary algebraic 
manipulations [ [Olmo et al., 200"9] ]. To do it, one defines a rank- two symmetric 
tensor W" v such that ^gg px JJr5\ + 2f Q B x u ) = y^hh 13 ", which turns (fl~35j) 
into the well-known equation [V— ~hh = 0, and implies that V^ v is given 
by the Christoffel symbols of the tensor h^ v , i.e., 

r /3 7 = — {dphp-t + d^hpp - dphpj) . (1.36) 

From the defining expression of h^, one finds that the relation between 
and g^v can be expressed as follows 

V = VdetlJp- 1 ]/^ , h» v = 9 ^£= , (1.37) 

VdetS 

where we have defined the matrix E a " = (/rS^ + 2fQB a u ). With these re- 
lations and definitions, the field equations for the metric can be written 
in compact form expressing (|1.34|) as B/En" = ^6" + K 2 Tp V and using the 

relation B/E/ = \/detT,R fla (h)h al/ to obtain [ |01mo fc Rubiera, 201 1] ] 

R/(h) = ^L=({-S; + ^TA . (1.38) 
VdetE \ 2 J 

In general, it will be more convenient to work with the field equations for the 
auxiliary metric h^ v because their form is more tractable. Nonetheless, if one 
insists on writing the field equations using the metric g^ v , one must note that 
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the connection (|1.36[) is related to the Levi-Civita connection of g^ by the 
tensor [recall Eq. (fl~T2|) ] 

A a M ee T$ 7 -L% = ^ [v£ h pv + V - V£ V] • C 1 - 39 ) 
The Riemann tensors of and Ui are thus related as follows 

i?V„(r) = R a j3^ v (L) + - + ^A"a - Kp A ^ > (i- 40 ) 

which allows to express (|1.38[) in terms of the Ricci tensor of the metric g^ v , the 
usual covariant derivatives of Lp , and the matter. 

1.4 f(R,Q) theories with a perfect fluid 

The explicit form of the matrix E that relates the metrics and <? M „ can 
only be found once all the sources that make up T^ u have been specified. In 
our discussion we will just consider a perfect fluid or a sum of non-interacting 
perfect fluids such that 

T„ v = (p + P)u„u v + Pg^ (1.41) 

with p — Pi and P = Pi. In order to find an expression for £, we first 
rewrite (|1.34[) using matrix notation as 

VqP 1 + SrB - f -I = n 2 f . (1.42) 
Using (|1.41[) this equation can be rewritten as follows 



Denoting A 2 = (k 2 P + j + an d making explicit the matrix representation, 
(|1.43p becomes 



2 



where /3X3 denotes 3-dimensional identity matrix. Since the right-hand side of 
(jl.44p is a diagonal matrix, it is immediate to compute its square root, which 
leads to 



where s± denotes a sign, which can be positive or negative, and 6*3x3 denotes 
a 3X3 diagonal matrix with elements {si = ±1}. For consistency of the theory 
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in the limit /g — > 0, we must have si = 1 and S3X3 — Isxs- This result allows 
to express £ as follows 

' fTi 

02/3x3 

where <7\ and (T2 take the form 



£ = 



(1.46) 



01 

CT 2 




2/ QA /A 2 - K 2(p + P) 
2fcA . 



(1.47) 



Note that we have kept the two signs ± in <j\ . The reason for this will be under- 
stood later, when particular models are considered. The point is that in some 
cases of physical interest, at high densities one should take the negative sign in 
front of the square root to guarantee that a\ is continuous and differentiable 
accross the point where the square root vanishes. This technical issue does not 
arise for 02- 

1.4.1 Workable models: f(R, Q) = f(R) + aQ 

So far we have made progress without specifying the form of the Lagrangian 
f(R, Q). However, in order to find the explicit dependence of R = B^ and Q — 
B^ a B a ^ with the p and P of the fluids, we must choose a Lagrangian explicitly. 
Restricting the function f(R, Q) to the family f(R, Q) = f(R) + aQ, we will see 
that it is possible to find the generic dependence of Q with p and P, while R 
is found to depend only on the combination T = —p + 3P [|01mo et al., 2009 ]. 
The reason for this follows from the trace of (|1.33|) with g^ v , which for this 
family of Lagrangians gives the algebraic relation RJr — 2/ = n 2 T and implies 
that R = R(T) (like in Palatini f(R) theories). For these theories, we have that 
fq = a, which is a constant. Therefore, from the trace of (|1.44l) we find 



(1.48) 



which can be cast as 




After a bit of algebra we find that 



A = 



3A 



X A -K 2 (p+P) 




(1.49) 



(1.50) 
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From this expression and the definition of A 2 , we find 




/q y /q 

(1.51) 



where R, /, and are functions of T = —p + 3P. 

1.5 Nonsingular cosmologies in f(R,Q) theories 



The difficulties faced by GR to provide a consistent description of singulari- 
ties and quantum phenomena at high energies (microscopic or Planck scales) is 
generally seen as an indication that we should go beyond the standard geomet- 
ric structures to successfully quantize the theory and avoid singularities. This 
idea has motivated a variety of approaches that range from the consideration of 
higher-dimensional superstrings and other extended objects [ |Green et al., 1987] ], 

to non-commutative geometries or non-perturbative quantization methods [| Ashtckar fc Lewandowski, 2004J ], 
[ |Rovelli, 2004] ] , [ [Thiemann, 2007] ] , to name just a few well-known cases. Un- 
fortunately, the formidable task of building a satisfactory quantum theory of 
gravity is not yet complete. Moreover, even if we managed to get such a theory, 
we would still have to face the challenge of testing its predictions. In this sense, 
it should be noted that since the quantum gravitational regime is so far from our 
current and future experimental capabilities, our only hope might be to use the 
information available in the cosmic microwave background radiation to verify 
or rule out our theories [ |Agullo fc Parker, 2011 ]. How much of the quantum 



gravitational regime could be contrasted with these yet-to-come theories is not 
clear. This is due, in part, because the theorized rapid accelerated expansion 
that took place during the inflationary period may have washed out many of 
the relevant proper signatures needed to distinguish the predictions of different 
quantum theories of gravity. 

A conservative approach, therefore, consists on exploring the quantum prop- 
erties and interactions of the matter fields in the very early universe using 
the well-established methods of quantum field theory in curved space-times 
[ |Parker &: Toms, 2009| ], The success of this approach has been confirmed in 
combination with models of inflation and sheds relevant light on the mechanisms 
that may have caused the primordial spectra of scalar and tensorial perturba- 
t ions [|Weinberg, 2008| ], [ |Lyth fc Liddle, 2008| ], [ |Dodelson, 2003| ],[ |Liddle fc Lyth, 2000| ], 
[ |Kolb fc Turner, T9 90 ]. The applicability of this approach, however, becomes 
unreliable at increasing energies as the regime of the classical big bang singu- 
larity is approached and the quantum fluctuations of the gravitational field can 
no longer be neglected. At that stage, a complete quantum theory of gravity 
seems necessary to provide a consistent description of the ongoing physical pro- 
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cesses. Obviously, different quantum theories could lead to completely different 
quantum gravitational scenarios and, therefore, a generic quantum origin for 
the universe cannot be guessed a priori by any logical means. 

In recent years, bouncing cosmological models have attracted much attention 
[ [Novello &: Per ez Bcrgliaffa, 2008']. These are scenarios in which the big bang 
singularity is replaced by a quantum-induced bounce that connects an earlier 
phase of contraction with the subsequent expanding phase (in which we happen 
to exist). In such scenarios, aside from the quantum regime, the contracting 
and expanding phases are expected to asymptote an effective classical geome- 
try whose dynamics, on consistency grounds, should match that of GR at low 
energies. In this context, and as an intermediate step between the quantum 
field theory approach in the (singular) curved background provided by GR and 
a (nonsingular) full theory of quantum gravity, one could consider the case of 
a smooth effective geometry free from big bang singularities on top of which 
quantum matter fields could still be treated perturbatively in a consistent way. 
This view would somehow disentangle the non-perturbative part of the quan- 
tum gravitational sector into an effective classical, nonsingular geometry, plus 
perturbative quantum corrections that propagate on top of the regular effective 
background. The absence of curvature singularities would make the treatment 
of quantum fields on the resulting geometry more reliable, and could help shed 
new light on the effects of the matter-gravity interaction in the very-early uni- 
verse. 



In the literature there exist many interesting examples of (quantum and non- 
quantum) cosmological models that avoid the big bang singularity by means of a 
bounce. Roughly, those models can be classified in two large groups, depending 
on whether they contain a modified gravitational sector or a modified matter 



sector (see [Novello fc Perez Bergliaffa, 2008 for details and a very complete 



list of references). Generically, modified gravity theories imply the existence 
of new dynamical degrees of freedom, such as gravitational scalar fields (like 
in scalar-tensor theories), higher-derivatives of the metric, extra dimensions, 
. . . The consideration of exotic matter sources may be justified, in some cases, 
from an effective field theory approach, such as in the case of non-linear theories 
of electrodynamics, which naturally arise in low-energy limits of string theories. 
In the remainder of this chapter, we are going to study bouncing cosmological 
models from the modified gravity perspective provided by the Palatini theories 
discussed above. This approach is particularly interesting because, despite be- 
ing a modified-gravity approach, the underlying mechanisms that modify the 
gravitational dynamics are not associated with new dynamical degrees of free- 
dom or higher-derivative equations. In fact, it is the nontrivial role played by 
the matter in the determination of the space-time connection that induces non- 
linearities in the matter sector that end up changing the dynamics at very high 
matter-energy densities. In this sense, it should be noted that the gravitational 
field equations in vacuum exactly recover those of GR (with possibly a cosmolog- 
ical constant, depending on the particular Lagrangian chosen). For this reason, 
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this type of theories can be regarded as a minimal extension of the standard 
model of gravitational physics, because they only appreciably depart from GR 
in regions that contain sources and when those sources reach the energy-density 
scales that characterize the correcting terms of the Lagrangian. 



1.5.1 Homogeneous cosmologies in f(R,Q) theories 

In this section we introduce the basic definitions and formulas needed to derive 
the equations for the evolution of the expansion and shear [[Wald, 1984 ] for 
an arbitrary Palatini f(R,Q) theory of the kind presented in Section 14.3.21 . 
These magnitudes will be very useful to extract information about the geomet- 
ric properties of the space-time and to determine whether cosmic singularities 
are present or not. We focus on homogeneous cosmologies of the Bianchi I 
type (a different expansion factor for each spatial direction) because that will 
allow us to test the rebustness of our results against deviations from the ideal- 
ized Friedmann-Robertson- Walker space-times (same expansion rate in all the 
spatial directions). We will also particularize our results to the case of f(R) 
theories, i.e., no dependence on Q. 



We consider a Bianchi I spacetime with physical line element of the form 

3 

ds 2 = g^dx^dx" = -dt 2 + ^ a 2 (t){dx 1 ) 2 (1.52) 



In terms of this line element, using the relation between metrics (|1.37l) and the 
expression (|1.46[) for the matrix £ of a collection of perfect fluids, the nonzero 
components of the auxiliary metric h^ v become 



( °l 



-S 



hij = \/o-\(i2a 2 5ij 
The relevant Christoffel symbols associated with h^ v are the following: 



fla 2 Sij 



r* 
1 u 



r*. 



tj 



25 
fiof 
25 

£ 
2 



ft 2d, 
1 

n 2d, 
o + 



The nonzero components of the corresponding Ricci tensor are 

3Q 3fl f S fl\ 1 ( S 2tl 



Rttih) 



2fl An \ S 



2 1 s n 



(1.53) 
(1.54) 

(1.55) 
(1.56) 

(1.57) 



15 



2 S 



2H 



2Hi i J2 H k 



ft ' 

if 3D, S' 
2 l 7T 



m 3n 
2~n \ IT 




(1.59) 



where Hk = dfe/flfe- For completeness, we give an expression for the correspond- 
ing scalar curvature 



R(h) = ~ 



2£> 



it 











(i- 





(1.60) 

From the above formulas, one can readily find the corresponding ones in the 
isotropic, flat configuration by just replacing ifj — > H. For the spatially nonflat 
case, the Rtt(h) component is the same as in the fiat case. The Rij(h) compo- 
nent, however, picks up a new piece, 2K^/ij, where "fij represents the nonflat spa- 
tial metric of — a^ij. The Ricci scalar then becomes R(h) — > R K=0 (h) + ^^. 



305 
2flS 



1.5.2 Shear 



From the previous formulas and the field equation (|1.38l) . we find that the 



combination Ri 



Rj 3 (no summation over indices) leads to 



s 



Hi 



Hi, 



= o , 



(1.61) 



where we have defined ify = Hi — Hj. Note that the final equality Ri~Rj J =0, 
follows from the fact that the right hand sides of E4 1 and Rj J as given by (|1.38|) 
are equal. Expressing (|1.61j) in the form 



R 1 - R j = — f 



In Hij + ln( ai a 2 a 3 ) + In O 3 / 2 - In S 1/2 







(1.62) 



we see that it can be readily integrated regardless of the number and particular 
equations of state of the fluids involved. The result is 



" "(>:: (a ia2 a 3 ) - at V(t) ' 



(1.63) 



where the constants C*. 



-Cji satisfy the relation C%2 + C23 + C31 = 0, 



Vq represents a reference volume, and V(t) = Voa\a2a 3 represents the volume 
of the universe. It is worth noting that writing explicitly the three equations 
(jl.63l) and combining them in pairs, one can write the individual Hubble rates 
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as follows 



(Cl2 — C31) , 


( V 


3(Ti 


W(t) 


(C23 — C12) 


( V 


3(Tl 


KV(t) 


(C31 — C23) 


( V 


3(Tl 


KV(t) 



H 2 = 9 + K *\ -£r (1.64) 



where 9 is the expansion of a congruence of comoving observers and is defined 
as W = J^i^i- Using these relations, the shear a 2 — ~}2i(Hi — 9) 2 of the 
congruence takes the form 

2 _ (C12 +g 2 2 3 + C|i) / Vb \ 2 , , 

where we have used the relation (C12 + C23 + C31) 2 = 0. 
1.5.3 Expansion 

We now derive an equation for the evolution of the expansion with time and a 
relation between expansion and shear. From previous results, one finds that 

CM , -i £ Hi + I fa H.) 2 + S £ * + I (2) 2 (1.66) 
In terms of the expansion and shear, this equation becomes 

G » s3 ( e + Ji) 2 -f < L67 » 

From the field equation f| 1 . 38[) . we find that 

ou = i±^±m , (1 . 68) 

which in combination with (jl.67p yields 

For a set of non- interacting fluids with equations of state Wi = Pi/pi, we have 
that fl = Q(pi,Wi) and, therefore, fl = X^^wPii where fl Pi = dfl/dpi. Since 
for those fluids the conservation equation is pi — —39(1 + u>i)pi, we find that 
(l = —30 5^(1 + Ui)pifl Pi . With this result, (|1.69l) can be written as 

w 1+§Al v /±^>4, (1 , 0) 
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where we have defined 



A 1 = -£(l + ^%^. (1.71) 

i 

Note that in this last equation Wi = Wi(pi), i.e., they need not be constants. 
For fluids with constant Wi , the conservation equation implies that their density 

depends on the volume of the universe according to pi(t) = pi(t Q ) [yf^j 
This implies that once a particular Lagrangian is specified, the equations of 
state Pi — WiPi are given, and the anisotropy constants Cy are chosen, the 
right-hand side of Eqs. (|1.65|) and (jl.70j) can be parametrized in terms of V(t). 
This, in turn, allows us to parametrize the Hi functions of (11.641) in terms of 
V(t) as well. This will be very useful later for our discussion of particular models. 

In the isotropic case (a 2 =0,6 = a/a = T-L) with nonzero spatial curvature, 
ljl.70jl takes the following form: 

H = ^ \ a - 1 1.72 

6^1 [l + |Ai] 2 

The evolution equation for the expansion can be obtained by noting that the 
Rij equations, which are of the form Rij = (Cl/2S)gij [•■•] = (/ /2 + n 2 P)gij /<T2, 
can be summed up to give 

1.5.4 Limit to f(R) 

We now consider the limit /q — > 0, namely, the case in which the Lagrangian 
only depends on the Ricci scalar R. Doing this we will obtain the corresponding 
equations for shear and expansion in the f(R) case without the need for extra 
work. From the definitions of A 2 (see below eq.([TT43|)), and o\ and 02 in (|1.47[) . 
it is easy to see that in the limit /q — > we get 

£71 -> cr 2 -> Ir (1-74) 
S -> -> f R . (1.75) 

With these rules it is easy to see that h^ u = fng^ , which makes (|1.38l) boil down 
to the expected field equations for Palatini f(R) theories, namely, fFt.R^ v {h) — 
^Qiiv = n 2 T^ u . Equation (| 1 .63[) turns into 
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from which one can easily obtain expressions for Hi, H2 and H3 as in fll .64[) . 
The shear becomes 



(C12 + C23 + C31) ( v o 



9/i 



(1.77) 



where C12 + C23 + C31 = 0. The relation between expansion and shear for a 
collection of non-interacting perfect fluids now becomes 



(1.78) 



where Ai is given by (jl.71l) but with SI replaced by Jr. In the isotropic case 
with nonzero K we find 



H 2 



1 



f + K 2 (p + 3P) 



GKf R 



QfE 



l + |Ai 



(1.79) 



1.5.5 Bouncing f(R) Cosmologies 

We now present the cosmological dynamics of simple f(R) models to illustrate 
how this family of theories modifies the standard Big Bang picture of the early 
universe. Consider, for instance, the modeH f(R) = R + aR 2 /R P , where R P = 
lp 2 = c 3 /HG is the Planck curvature. From the trace equation Rfp — 2/ = k 2 T 
(see Sec ll.4.1]) . we find that this model leads to the same relation between the 
matter and the scalar curvature as in GR, namely, R = —k 2 T. This implies 
that the theory behaves as GR whenever the energy density is much smaller 
than the Planck density scale pp = Rp/n 2 . Since by definition 6 — -| J2i = 

Inaia2a3 = ^y, where V = Voaia2Ct3 represents the volume of the universe 
(with Vo = V(to)), Eq. (|1.78p for this quadratic model with dust and radiation 
leads to 



Pd + Pr + 



2p P 



1 



2apd 
Pp 



(C 12 + C| 3 + Cfj) f V \ 



3(1-2^ 
' pp 



54 ( 1 - y& 
' pp 



2 {V(t)j ' 
(1.80) 

where p d = p dfl {yft^j and Pr = Pr,o 

In general, an homogeneous cosmological model experiences a bounce when the 
expansion 8 vanishes, which implies an extremum (a maximum or a minimum) 
of the volume of the Universe. If V(t) vanishes at some finite time, then a 
big bang or big crunch singularity is found, depending on whether V > or 
V < at that time. Focusing for the moment on the isotropic case, Cf 2 + 



2 Note that the constant a could be absorbed into a redefinition of Rp and, therefore, only 
its sign is relevant. 
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Figure 1.1: Representation of the Hubble function (left) and volume of the 
Universe (right) as a function of time for the model f(R) = R — R 2 /2Rp in 
a universe filled with dust and radiation (for the numerical integration pdfi = 
I0 3 p r ,o, and V = 10 5 Vb). The GR solutions corresponding to a contracting 
branch, which ends in a big crunch, and an expanding branch, which begins 
with a big bang, are represented together with the bouncing solution of the 
Palatini model that interpolates between those singular solutions. 



C| 3 + Cfj = 0, we find that a bounce occurs if a < when pd reaches the value 
pf = pp/(2\a\) [see Fig ll.lj) ] . This value of the density implies that /r = 1 — 
2apd/pp = 0. This condition, fp = 0, characterizes the location of the bounce 
in Palatini f(R) theories with a single fluid with constant equation of state 
[ |Barragan fc Olmo, 2010] ]. For our quadratic model, in particular, bouncing 
solutions exist if the dynamics allows to reach the density ps = 2a(3w-i) > 0- 
This means that for a > fluids with w > 1/3 avoid the initial singularity, 
whereas for a < it takes w < 1/3. The case a = naturally recovers the 
equations of GR. It is worth noting that a cosmic bounce may arise even for 
presureless matter, w — 0, if a < 0, which implies that exotic sources of matter- 
energy that violate the energy conditions are not necessary to avoid the big 
bang singularity in this framework. The reason for this is that at high energies 
gravitation may become repulsive for matter sources with w > — 1, whereas it 
is attractive at low energy-densities for those same sources. 
Note also that the pure radiation universe, w = 1/3, is a peculiar case because it 
does not produce any modified dynamics in Palatini f(R) theories. On physical 
grounds, however, it should be noted that due to quantum effects related with 
the trace anomaly of the electromagnetic field, a gas of photons in a SU(N) 
gauge theory with Nf fermion flavors has an effective equation of state given by 



w 



rad 1 5a 2 (JV e + fJV / )(^JV e -|JV / ) 
eft 3 18^ 2 + 



(1.81) 



where N c is the color number of the gauge theory (which has N C (N C — 1) gener- 
ators) [ jKajante et al., 2003| ], [ |Lambiase fc Mohanty, 2007] ] . Therefore, a uni- 
verse filled with photons should be able to avoid the singularity if a > 0. 
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Anisotropic Bounce in f(R)=R - jf- Anisotropic Bounce in f(R)=R - £ 




Figure 1.2: Representation of the expansion (left) and volume of the Universe 
(right) as a function of time for the model f(R) — R—R 2 /2Rp in a universe filled 
with dust and radiation with anisotropies (for the numerical integration pdfi = 
10 3 Pr,0: an d V = 10 5 Vo). From right to left, we have plotted the bouncing cases 
C 2 = 0, 40, 40.60211073, 40.60211073942454489657, and the collapsing case with 
C 2 = 40.60211073942454489658. Fine tunning the value of C 2 even more should 
allow to keep the universe in its minimum for longer periods of time in the past, 
which eventually should lead to an asymptotically static solution. 

In physically realistic scenarios, one should consider the co-existence of sev- 
eral fluids and take into account the time dependence of the number of ef- 
fective degrees of freedom and the transfer of energy among different species 
[ |Kolb fc Turner, 1990| 1, which leads to the possibility of having different effec- 
tive fluids at different stages of the cosmic expansion. In this sense, the "dust 
plus radiation" model represented by (|1.80[) needs not be accurate at all times 
because dust particles may become relativistic at high energies and contribute 
to p r rather than to pd- This suggests that the choice/determination of the sign 
of the parameter a is not a trivial issue and would require a very careful and 
elaborate analysis (which goes beyond the scope of this introductory work). 

When anisotropies are taken into account, one finds that bouncing solutions 
are still possible as long as the amount of anisotropy is not too large. In Fig |1.2[ 
we see that increasing the value of C 2 = C\ 2 + C 2 3 + C 31 from zero, the vol- 
ume of the universe presents a minimum as long as C 2 < C 2 . If C 2 > C 2 , the 
collapse is unavoidable and V — > in a finite time. The critical case C 2 — > C 2 
represents a configuration that is neither a bouncing universe nor a big bang. It 
corresponds to a state in which the volume of the universe remains constant in 
the past and expands in the future. Though this solution is clearly unstable and 
fine-tuned, its existence puts forward the possibility of obtaining static regular 
solutions corresponding to ultracompact objects, which could shed new light on 
the internal structure of black holes and/or topological deflects when Planck 
scale corrections to the gravitational action are taken into account. It should be 
noted, however, that in order to obtain this asymptotically static solution one 
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must cross from the domain where fp > to the region where fp<0. Since the 
shear, as defined in (jl.77[) for f(R) theories with perfect fluids, is proportional 
to 1//^, the crossing through — implies a divergence in some curvature 
scalars of the theory. Whether this divergence is a true (or strong) physical 
singularity in the sense defined in [ |Tipler, 1977| ], [ |Clarke &: Krolak7 l985 ], 
[ |Krolak,T"9 98 ] is an open question that will be explored elsewhere. In any 
case, we remark that the existence of that divergence does not have any effect 
on the time evolution of the expansion 9, as can be seen in Fig |1.2l 



1.5.6 Nonsingular Universes in f(R,Q) Palatini theories 

In the previous section we have seen that Palatini f(R) models are able to avoid 
the big bang singularity in idealized homogeneous and isotropic scenarios but 
run into trouble when anisotropies are present. The divergence of the shear is 
a generic problem for those f(R) theories in which the function fp vanishes 
at some point, regardless of the number and equation of state of the fluids 
involved. Though the nature of this divergence has not been identified yet with 
that of a strong singularity, which besides the divergence of some components 
of the Riemann, Ricci, and Weyl tensors also requires the divergence of some 
of their integrals, its very presence is a disturbing aspect that one would like to 
overcome within the framework of Palatini theories. In this sense, a natural step 
is to study the behavior in anisotropic scenarios of some simple generalization 
of the f(R) family to see if the situation improves. Using Lagrangians of the 
form presented in (|1.4.1[) . we will show next that completely regular bouncing 
solutions exist for both isotropic and anisotropic homogeneous cosmologies. 



Isotropic Universe 

Consider Eq. (|1.72j) particularized to the following f(R,Q) Lagrangian 

f(R,Q) = R+a^- + b-^- (1.82) 
rip rip 

For this theory, we find that R = n 2 (p — 3P) and Q = Q(p,P) is given by 
(ll.51|) with a = b/Rp. From now on we assume that the parameter b of the 
Lagrangian is positive and has been absorbed into a redefinition of Rp, which 
is assumed positive. This restriction is necessary (though not sufficient) if one 
wants the scalar Q to be bounded from above when fluids with w > — 1 are 
considered. Stated differently, when b/Rp > 0, positivity of the square root of 
Eq. (|1.5ip establishes that there may exist a maximum for the combination p+P. 

In order to have (|1.72l) well defined, one must make sure that the choice of 
sign in front of the square root of o\ in (ll.47|) is the correct one. In this sense, 
we find that to recover the f(R) limit and GR at low curvatures, we must take 
the positive sign, i.e., o\ = o\ ■ However, when considering particular models, 
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Discontinuity of o~i + and o~j 

ft 




0.0 " 0.1 0.2 0.3 0.4 P 



Figure 1.3: Illustration of the need to combine the 
two branches of a\ to obtain a continuous and dif- 
ferentiable curve. The branch that starts at <j\ — 1 
has the plus sign in front of the square root (con- 
tinuous green line). When the square root vanishes 
(at the blue dot), the function must be continued 
through the dashed red branch, which corresponds 
to the negative sign in front of the square root. 



which are characterized by the constant a and, for instance, a constant equa- 
tion of state w, one realizes that the square root may reach a zero at some high 
density. Beyond that point, we may need to switch from <r^ to to guarantee 
that <7i is a continuous and differentiable function (see Fig l 1.31 for an illustration 
of this point). Bearing in mind this technical subtlety, one can then proceed 
to represent the Hubble function for different choices of parameters and fluid 
combinations to determine whether bouncing solutions exist or not. 

The classification of the bouncing solutions of the model (|1.82p with a fluid 
with constant w was carried out in [[Barragan & Olmo, 2010 ] . It was found that 
for every value of the parameter a there exist an infinite number of bouncing 
solutions, which depend on the particular equation of state w. The bouncing 
solutions can be divided into two large classes: 

• Class I: a > 0. The bounce occurs when the scalar Q reaches its maxi- 
mum value and happens for all equations of state satisfying the condition 

a 

W > Wmin = - ; - • (1.83) 
2 + 3a 

From this equation it follows that a radiation dominated universe, with 
w = 1/3, always bounces for any a > 0. 
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• Class II: a < 0. This case is more involved because the bounce can occur 
either at the point where Q reaches its maximum or when a\ vanishes. 
This last case can only happen at high curvatures when we are in the 
branch defined by o\ = rxf. To proceed with the classification, we divide 
this sector into several intervals: 

— If —1/4 < a < 0. The bounce occurs if 



We see that when a = we find agreement with the discussion of case 
I. As a approaches the limiting value —1/4, the bouncing solutions 
extend up to w — > —1/3. 

— If— 1/3 < a < —1/4. Numerically one finds that the bouncing 
solutions cannot be extended below w < — 1 and occur if — 1 < w < 

00, where w = — 1 is excluded. 

— If— l<a<— 1/3. In this case, one finds numerically that the 
bouncing solutions are restricted to the interval — 1 < w < r^tti > 

1, where a = 1.1335 and f3 = -3.3608. 

— If a < — 1. Similarly as the family a > 0, this set of models also 
allows for a simple characterization of the bouncing solutions, which 
correspond to the interval — 1 < w < a/(2 + 3a). In the limiting case 
a = — 1 we obtain the condition — 1 < w < 1 (compare this with the 
numerical fit above, which gives — 1 < w < 1.12). 

1.5.7 Anisotropic Universe 

Using Eqs. (|1.70[) and (|1 .72|) . the expansion can be written as follows: 



where H represents the Hubble function in the K = isotropic case. To bet- 
ter understand the behavior of (9 2 , let us consider when and why H 2 vanishes. 
Using the results of [ |Barragan fc Olmo, 2010| ] summarized above, one finds 
that H 2 vanishes either when the density reaches the value PQ max or when the 
function ct\ vanishes. These two conditions imply a divergence in the quantity 
(1 + |Ai) 2 , which appears in the denominator of H 2 and, therefore, force the 
vanishing of H 2 (isotropic bounce). Technically, these two types of divergences 
can be easily characterized. From the definition of Aj in (11.711) . one can see that 
A 1 ~ dpCl/Cl. Since 51 = Ja\a2, it is clear that Ai diverges when a\ = 0. The 
divergence due to reaching PQ mox is a bit more elaborate. One must note that 
dpil contain terms that are finite plus a term of the form d p X, with A defined 
below Eq. (|1.43[) . In this A there is a Q term hidden in the function f(R,Q), 
which implies that d p X ~ d p Q/Rp plus other finite terms. From the definition 




(1.84) 
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Figure 1.4: Representation of the expansion squared (left) and volume of the 
Universe (right) as a function of time for the model f(R, Q) = R — lpR 2 /2 + lpQ 
in a radiation universe with anisotropies. We have plotted the bouncing cases 
C 2 = 0, 50, 10 2 , 10 3 . Note that the bounce always occurs at the same maximum 
density (minimum volume) . Note also that the time spent in the bouncing region 
decreases as the anisotropy grows. The starting point of the time integration is 
chosen such that at t = the two branches of erj coincide. 

of Q it follows that d p Q has finite contributions plus the term i9 p $/\/$, where 
$ = (1 + (1 + 2a)R/R P ) 2 - 4k 2 (p + P)/R P , which diverges when $ vanishes. 
This divergence of d p Q indicates that Q cannot be extended beyond the maxi- 
mum value Qmax- 

Now, since the shear goes like er 2 ~ 1/(cti) 2 [see Eg. (11.651) ]. we see that the 
condition o\ — implies a divergence on a 2 (though 2 remains finite). This is 
exactly the same type of divergence that we already found in the f(R) models, 
where <j\ — > fp. Since in the f{R) models the bounce can only occur when 
fit = 0, there is no way to avoid the divergence of the shear in the anisotropic 
case within the f(R) setting. On the contrary, since the quadratic f(R,Q) 
model (|1.82p allows for a second mechanism for the bounce, which takes place 
at PQ mal , there is a natural way out of the problem with the shear. 
Summarizing, we conclude that for universes governed by the Lagrangian (jl.82l) 
and containing a single stiff fluid there exist completely regular bouncing so- 
lutions in the anisotropic case for w > 2 +3a if a > 0, for wq < w < oo if 
-1/3 < a < 0, for w < w < (a + /3a)/(l + 3a) 2 if -1 < a < -1/3, and for 
— 1/3 < w < a I (2 + 3a) if a < — 1, where wo < is defined as the equation 
of state for which the (isotropic) bounce occurs when Q — Q max and a\ = 
simultaneously (see |Barragan fc Olmo, 2010] for details). These results imply 
that for a < the interval < w < 1/3 is always included in the family of 
completely regular isotropic and anisotropic bouncing solutions, which contain 
the dust and radiation cases. For a > 0, the radiation case is always nonsingular 
too. 
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1.5.8 Example: radiation universe. 



As an illustrative example, we consider here the particular case of a universe 
filled with radiation. Besides its obvious physical interest, this case leads to 
a number of algebraic simplifications that make more transparent the form of 
some basic definitions 



a 2 
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(1.86) 

(1.87) 

(1.88) 
(1.89) 



Note that the coincidence of the two branches of u\ occurs at n 2 p — Rp/6, 
where af — §. It is easy to see that at low densities (|1.86[) leads to Q « 



2 ■ 

4(«; 2 p) 2 /3 + 32(K 2 p) 3 /9i?p + 320(K 2 p) 4 /27R 2 p + . . ., which recovers the expected 
result for GR, namely, Q = 3P 2 + p 2 . From this formula we also see that the 
maximum value of Q occurs at K 2 p ma x — 3Rp/l6 and leads to Q ma x = 3R p /16. 
At this point the shear also takes its maximum allowed value, namely, cr^^ = 
yJiJlQRy, 2 (C 2 2 + C| 3 + Cf-J, which is always finite. At p max the expansion 
vanishes producing a cosmic bounce regardless of the amount of anisotropy [see 

Fignai. 



1.6 Conclusions and open questions 

In this chapter we have tried to convey the idea that in the construction of 
extended theories of gravity, one should bear in mind the fact that metric and 
connection are equally fundamental and independent objects. This observation 
allows to broaden the spectrum of available possibilities to go beyond the stan- 
dard model of gravitation. In fact, any theory of gravity based on a geometry in 
which the connection has been forced to be given by the Christoffel symbols of 
the metric admits an alternative formulation in which the form of the connec- 
tion is dictated by the theory itself, i.e., it is not given by convention or selected 
on practical grounds. 

In our exploration of Palatini theories, we have seen that assuming that met- 
ric and connection are independent geometrical objects has non-trivial effects on 
the resulting field equations as compared with the usual metric formulation of 
the same theories. For the particular family of f(R, Q) models studied here, we 
have seen that the metric is governed by second-order equations that boil down 
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Dust+Radiation -Vs- Pure Radiation 
C 2 =10 3 
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Figure 1.5: Comparison of the expansion in a 
universe filled with dust and radiation (po.rad = 
10~ 3 po,dust) and a radiation dominated universe 
(dashed lines) for several values of the anisotropy. 



to GR in vacuum. This is in sharp contrast with the usual metric formulation 
of those same theories, where one finds fourth-order derivatives of the metric 



(see, for instance, [ Anderson, 1983 ] for a detailed analysis of the cosmology of 
the quadratic model (|1.82[) in metric formalism). The absence of higher-order 
derivatives in the Palatini formulation is a remarkable point that seems not to 
have been sufficiently appreciated in the literature. In fact, having second-order 
field equations is very important because it automatically implies the absence of 
ghosts and other dynamical instabilities. In this sense, it should be noted that 
LoveloclH theories [ |Zanelli, 2005 ], which are generally regarded as the natural 



extension of the Einstcin-Hilbert Lagrangian to higher dimensions, have received 
a lot of attention in the literature because they are seen as the most general 
actions for gravity that give at most second-order field equations for the metric. 
As we have seen here, this property is shared (at least) by all Palatini theories 
of the f(R,Q) type (with or without torsion). This puts forward that Palatini 
theories, are natural candidates to explore new dynamics beyond GR. 

Before concluding, we would like to stress the fact that the quadratic Palatini 
model (|1.82p is able to avoid the big bang singularity in very natural situations, 
such as in pure radiation, pure dust, or dust plus radiation universes with or 
without anisotropics (see Fig ll.5[) . This observation has been possible thanks 
to the formulas presented in section II. 5[ where we have extended the analysis 
carried out in [ Barraga n fc Olmo, 2010| ] for a single perfect fluid with constant 



3 We would like to mention that when Lovelock theories are formulated a la Palatini, the 
resulting field equations are exactly the same as one finds in their usual metric formulation 



[ |Borunda et al., 2008] ]. 
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equation of state to include several perfect fluids with arbitrary equation of state 
w(p). This allows to explore the dynamics of realistic cosmological models with 
several fluids and is a necessary step prior to the consideration of the growth 
and evolution of inhomogeneities in these nonsingular backgrounds. 
Though the model (|1.82j) has been proposed on grounds of mathematical sim- 
plicity and motivated by the form of the effective action provided by per- 
turbative quantization schemes in curved backgrounds, its ability to success- 
fully deal with cosmological [ Barragan fc Olmo, 2010| ] and black hole singu- 



larities [ |01mo fc Rubiera, 201 lj ] as well as other aspects of quantum gravity 



phenomenology [ Olmo, 2011c ] demands further theoretical work to provide a 
more solid ground to it. In this sense, we note that the effective dynamics of 
loop quantum cosmology [j Ashtekar fc Singh, 201 lj ] in a Friedmann-Robertson- 
Walker background filled with a massless scalar can be exactly reproduced by 
a Palatini f(R) theory [|01 mo &: Singh, 2009 ;]. The extension of that result to 
more general space-times and matter sources could shed new light on the po- 
tential relation of (II .82)) with a more fundamental theory of quantum gravity. 
All these open questions will be considered in detail elsewhere. 

This work has been supported by the Spanish grants FIS2008-06078-C03-02, 
FIS2011-29813-C02-02, the Consolider Program CPAN (CSD2007-00042), and 
the JAE-doc program of the Spanish Research Council (CSIC). 



2N 



Bibliography 



[Agullo & Parker, 2011] Agullo, I. and Parker, L. , Gen. Rel. Grav. 43, 2541 
(2011) [Int. J. Mod. Phys. D 20, 2861 (2011)] [arXiv: 1106.4240 [astro- 



[Anderson, 1983] Anderson, P., Phys. Rev. D 28, 271 (1983). 

[Ashtekar & Lewandowski, 2004] Ashtekar, A. and Lewandowski, J. , Class. 
Quant. Grav. 21 (2004) R53. 

[Ashtekar & Singh, 2011] Ashtekar, A. and Singh, P., Class. Quant. Grav. 28, 
213001 (2011) [arXiv:1108.0893 [gr-qc]]. 

[Barragan & Olmo, 2010] Barragan, C. and Olmo, G. J. , Phys. Rev. D 82, 
084015 (2010), [arXiv:1005.4136 [gr-qc]]. 

[Birrell & Davies, 1982] Birrel N.D. and Davies P.C.W., Quantum fields in 
curved space, (Cambridge University Press, Cambridge, England, 1982). 

[Borunda et al., 2008] Borunda, M., Janssen, B. and Bastero-Gil, M., JCAP 
0811, 008 (2008) [arXiv:0804.4440 [hep-th]]. 

[Clarke & Krolak, 1985] Clarke, C. J. S. and Krolak, A. , Journ. Geom. Phys. 
2, 17 (1985). 

[Dodclson, 2003] Dodclson S., Modern Cosmology, Academic Press, (2003). 

[Feigl & Brodbeck, 1953] Feigl, H. and Brodbeck, M. (eds.), Readings in the 
phylosophy of science, Applcton-Century-Crofts, New York (1953). (See 
page 193). 

[Green et al., 1987] Green, M. , Schwarz, J. , and Witten, E. , "Superstring 
Theory", Cambridge University Press, Cambridge, England (1987). 

[Hchl et al., 1995] Hehl, F. W., McCrea, J. D., Mielke, E. W., and Ne'eman, 
Y., Phys. Rept. 258, 1 (1995) [gr-qc/9402012]. 

[Kajante et al., 2003] Kajantie, K., Laine, M. , Rummukaincn, K. and 
Schroder, Y. , Phys. Rev. D 67, 105008 (2003). 



ph.CO] 



29 



[Kibble, 1961] Kibble, T. W. B., J. Math. Phys. 2, 212 (1961). 

[Kolb & Turner, 1990] Kolb, E.W. and Turner, M.S. The early universe, West- 
view Press (1990). 

[Krolak, 1998] Krolak, A. , Class. Quant. Grav. 3, 267 (1998). 

[Lambiasc & Mohanty, 2007] Lambiase, G. , Mohanty, S. , JCAP 0712 (2007) 
008. 

[Liddlc & Lyth, 2000] Liddle A. R. and Lyth D.H., Cosmological inflation and 
large-scale structure, Cambridge University Press, (2000). 

[Lyth & Liddle, 2008] Lyth D.H. and Liddlc A.R. The primordial density per- 
turbation, Cambridge University Press, (2009) 

[Nakahara, 1990] Nakahara, M. , Geometry, Topology and Physics, Adam 
Hilger, Bristol (1990). 

[Novello & Perez Bergliaffa, 2008] Novcllo, M. and Perez Bcrgliaffa, S.E. , 
Phys.Rep. 463, 127-213 (2008). 

[Olmo & Singh, 2009] Olmo, G. J. , and Singh, P., JCAP 0901, 030 (2009) 
[arXiv:0806.2783 [gr-qc]]. 

[Olmo et al., 2009] Olmo, G. J. , Sanchis-Alepuz, H., and Tripathi, S., Phys. 
Rev. D80, 024013 (2009). [arXiv:0907.2787 [gr-qc]]. 

[Olmo, 2011a] Olmo, G. J. , Int. J. Mod. Phys. D 20, 413 (2011). 

[Olmo, 2011b] Olmo, G. J. , AIP Conf. Proc. 1458, (2011) 222-237. 

[Olmo, 2011c] Olmo, G. J. , JCAP 1110, 018 (2011) [arXiv: 1101. 2841 [gr-qc]]. 

[Olmo & Rubiera, 2011] Olmo, G. J. , and Rubiera-Garcia,D. , Phys. Rev. D86 
(2012) 044014; Int.J.Mod.Phys. D21 (2012) 1250067; Eur.Phys.J. C72 
(2012) 2098. 

[Olmo & Rubiera, 2012] Olmo, G. J. , and Rubiera-Garcia,D. , to appear 
(2012). 

[Parker & Toms, 2009] Parker L. and Toms D.J., Quantum field theory in 
curved spacetime: quantized fields and gravity, (Cambridge University 
Press, Cambridge, England, 2009). 

[Poplawski, 2010] Poplawski, N. J., Phys. Lett. B 694, 181 (2010) [Erratum- 
ibid. B 701, 672 (2011)] [arXiv: 1007.0587 [astro-ph.CO]]. 

[Rovclli, 2004] Rovelli, C. , Quantum Gravity, (Cambridge U. Press, 2004). 

[Schouten, 1951] Schoutcn, J. A. , Tensor analysis for physicists, Oxford Uni- 
versity Press, London (1951). 



30 



[Sciama, 1962] Sciama, D. W. , in: Recent Developments in General Relativity, 
p. 415 (Pergamon, 1962). 

[Sciama, 1964] Sciama, D. W., Rev. Mod. Phys. 36, 463 (1964). 

[Thiemann, 2007] Thiemann, T. .Modern canonical quantum general relativity, 
(Cambridge U. Press, 2007). 

[Tipler, 1977] Tipler, F. J. , Phys. Lett. A 64, 8 ( 1977). 

[Wald, 1984] Wald, R.M., General Relativity, (The University of Chicago Press, 
Chicago, 1984). 

[Weinberg, 2008] Weinberg S., Cosmology, Oxford University Press, (2008). 

[Will, 1993] Will, CM. (1993). Theory and Experiment in Gravitational 
Physics, Cambridge University Press, Cambridge; Living Rev.Rel. 
9,3,(2005), gr-qc/0510072 . 

[Zanelli, 2005] Zanelli, J. (2005). arXiv:hep-th/0502193v4 . 



31 



